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The new analysis variable £, shown to be powerful on the data taken with the final 
configuration of CELESTE, has been applied to data taken with previous detector 
configurations. First, the analysis is validated on Crab observations, and then the 
cuts for the blazar Mrk 501 are optimized using Mrk 421 data since the sources have 
similar declinations. Data from Mrk 501 was recorded in 2000 and 2001. The old 
analysis gave a 2.5 a excess. We obtain an excess of 2.9 a during this time and of 
4.9 a during May and June 2000 that we interpret as a 7-ray signal from Mrk 501, 
for which we calculate a flux of (6.9 ± 2.2) x 10 -7 photons m -2 s _1 . An upper 
limit from the other data with no signal is determined. 



1 Introduction 

CELESTE (Cherenkov Low Energy Sampling and Timing Experiment) was a Che- 
renkov experiment using the heliostats of the former Electricite de France solar plant 
in the French Pyrenees at the Themis site. It detected Cherenkov light from showers 
produced in the atmosphere by cosmic rays and the 7-rays coming from high energy 
astrophysical sources. The light is reflected to secondary optics and photomultipliers 
installed at the top of the tower. Finally it is sampled to be analysed 1 5 1. 

Two states of the experiment have to be distinguished to classify CELESTE data. 
During the first one (between September 1999 and June 2001), 40 heliostats were used 
with two types of pointing (single pointing: all heliostats at 11 km, double pointing: 
half at 11 km, half at 25 km), and during the second one (between September 2001 
and June 2004), 53 heliostats pointing at 11 km were used (of which 12 veto heliostats 
aimed wide for proton rejection). An analysis improvement was made using 53 he- 
liostat data for the Crab Nebula, a bright and stable source, the new analysis variable 
provided a sensitivity of 5.7 a/Vh on 4.7 h data, whereas the old data analysis gave 
2.0 and 3.4 a/Vh for single and double pointing [4|. The new analysis has been 
tested on 40 heliostat Crab data and gives better sensitivity as will be shown. 

CELESTE has taken 40 heliostat data on the blazar Mrk 501. The old analysis 
gave 2.5 a |2|. We present here the results obtained with the new analysis after opti- 
mizing cuts on the blazar Mrk 42 1 which has nearly the same declination as Mrk 501. 
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2 New Analysis applied to 40 Heliostat Crab Data 



After data selection based on stability criteria, a software trigger is applied to the 
data in order to remove trigger bias. The data can then be analysed with the new 
variable £ which is a normalised sum of the Cherenkov pulses from each heliostat, 
optimised by finding the shower position that gives the narrowest sum (i.e. the most 
gamma-like) |3][Q. £ has been developed on Crab data taken with 53 heliostats. 
The 40 heliostats data sets have been analysed with new cuts on £, determined from 
simulations and data. The study confirmed the dependence of £ on configuration 
(heliostat number, trigger) and pointing (hour angle). Different cuts were chosen 
depending on these parameters, yielding a sensitivity of 4.0 u/yh and 5.2 er/Vh for 
single and double pointing with 40 heliostats. 

The acceptances of the detector were calculated to determine a flux. They depend 
on many parameters: atmosphere, optical simulation (heliostats, secondary mirrors), 
electronic simulation, and pointing (configuration, declination, hour angle). The si- 
mulations have been improved to better fit the data, taking into account variations of 
the trigger rates between different data sets (two different atmospheres) and degra- 
dation of the optical chain of the detector. The lightcurve for the Crab nebula was 
determined for the 40 heliostat data with a spectral hypothesis of E a+I3 log E with E 
in TeV, a = —2.74 and (3 — —0.50 |4|. It is shown in figure^ where the calculated 
flux is stable, within the statistical error bars. 
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Figure 1: Crab nebula lightcurve seen by CELESTE above 100 GcV with statistical error bars (monthly 
averages centred on New Moon). The dashed line represents the mean integral flux for all the data: 
4>(E > 100 GeV) = (3.6 ± 0.2) X 10~ 6 photons m~ 2 s" 1 ). MJD (Modified Julian Date) 51450 
and 51900 correspond to 09/29/1999 and 12/22/2000. 
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3 Mrk 501 Flux Determination after Cut Optimisation 
using Mrk 421 Observations 



The blazar Mrk 421 has the double advantage of having about the same declination 
as Mrk 501 and of being well observed with CELESTE. So cuts were optimized on 
Mrk 421 data. (In fact, both the simulation and the data lead to very nearby the same 
cuts.) We obtained a final significance for all Mrk 421 data of 29.4 a during 38.9 h. 
The integral flux is 4>(E > 100 GeV) = (2.3 ± 0.1) x 10~ 6 photons ixr 2 s" 1 for an 
E~ 2 spectral hypothesis. 

The number of events and significance for Mrk 501 after the same cuts are shown 
in tabled A 2.9 a excess is obtained for the 2000-2001 period. The lightcurve of 
Mrk 501 (figure |2 left top), determined for an E~ 2 spectral hypothesis, presents an 
excess of activity during May and June 2000 (between MJD 51660 and 51700). For 
this period, the significance of the excess is 4.9 a which is interpreted as a 7-ray signal 
from Mrk 501, since the blazar is already well-known in low- and high-energy 7-rays. 
Thus, a differential flux is calculated for this period and an upper limit is determined 
for the rest of the data which contains no signal. 



Cut 


Number of events 


Significance 


Signal to noise ratio 




N 0N 


NoFF 


NoN — NoFF 


N a 


NoN — NoFF rnyl 
N [%1 

NOFF 


Raw trigger 


815 924 


813 641 


2 284 


1.5 


0.3 


Software trigger 


472450 


469 569 


2 881 


2.5 


0.6 


Final cuts 


10 565 


10 096 


469 


2.9 


4.6 



Table 1: Number of events for the blazar Mrk 501 after cuts for the final data set (11.5 h). 



4 Discussion 

The differential flux and the upper limit, obtained for a counting maximum at 1 10 GeV, 
are drawn on figure[2]right. The different states of activity of Mrk 501 in 7-rays are not 
correlated to variations in X-rays: in figure[2]left bottom, the ASM flux is 40 ± 10 % 
higher in May and June 2000 than the low baseline of the rest of the studied period 
(figure [2] right shows a x20 and xlOO increases during historical flares. These ob- 
servations suggest that for May and June 2000, the emission population was different 
between synchrotron emission in X-rays and 7-ray emission. If the 7-ray emission is 
due to an inverse Compton process, it can not be explained with a SSC (Synchrotron 
Self-Compton) model for this particular case. Further observations of Mrk 501 with 
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Figure 2: Left: Mrk 501 lightcurves (monthly averages centred on New Moon) seen by CELESTE 
Q>(E > 110 GeV) = (6.9 ±2.2) X 10~ 7 photons m~ 2 s"\ top) and RXTE/ASM in X-rays (f = 
0.40 ± 0.01 photons s" 1 , bottom) with statistical error bars. MJD 51600 and 52000 correspond to 
02/26/2000 and 04/01/2001. Right: overall spectral energy distribution of Mrk 501 of 04/16/1997, 
04/29/1998 and June 1999 |6|. Red circle (with statistical and systematic error bars) and arrow are the 
CELESTE flux and upper limit determined in this work at 2.6 X 10 25 Hz. 

new generation detectors (GLAST, HESS) and simultaneous to other wavelength ob- 
servations will help to understand the emission processes of Mrk 501. 
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